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The current status of the application of hydrodynamics to ultrarelativistic heavy 
*•) ion collisions is reviewed. We elaborate on the arguments for strong transverse 

^j- _ flow and rapid thermalization and discuss future applications and trends in hydro- 

f*^ , dynamics. 

1. Overview 

Q i 

^ | Among the many predictions 1 made for obscrvables of collisions at the Rel- 

Ph ■ ativistic Heavy Ion Collider (RHIC), hydrodynamic calculations regarding 

the expected flow characteristics 2,3 turned out to be surprisingly accurate 
already a few months after running the new facility. 4 In fact, those predic- 
tions were so convincing that they were awarded a bottle of wine in the 
RHIC prediction competition. 5 That RHIC creates systems which - for the 
dominant part of their lifetimes - expand according to hydrodynamic prin- 
ciples is of utmost importance: The dynamical evolution is governed by 
the nuclear equation of state at extreme energies, whose thermodynamic 
features are thus accessible through collisions in the laboratory. 

In the following, I will review indications of the strong transverse ex- 
pansion that can be deduced from the observable final state of central and 
non-central Au+Au collisions. Thereby, I will pay particular attention to 
signals of collectivity at mid-rapidity, i.e. in the plane perpendicular to the 
beam-axis in the center of mass of the reaction (at vanishing pseudorapid- 
ity rj). Although interesting phenomena are to be expected at 77 ^ 0, to 
which I will return later, the midrapidity region is of particular interest for 
the observation of collectivity: Before the reaction, there is no transverse 
component of motion. Observing collectivity in the transverse plane at 
midrapidity must thus be a result of the dynamical evolution of the fireball 
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after the initial impact. In contrast, at rapidities r\ ^ 0, it is difficult to 
disentangle collectivity created by the fireball region and the longitudinal 
dynamics of the nuclei that was present already before the reaction. 

To address the questions of the how and when those collective features 
in the transverse plane are established, I will turn to dynamical models 
and specifically introduce the formalism of hydrodynamics to study the 
temporal evolution of the fireball matter in this framework. This analysis 
makes the success of the simple formulas of blast wave parameterizations 
understandable and validates such a simplistic approach as a tool to rapidly 
characterize and summarize the main thermodynamic information of the 
final state. 

I will then highlight the most important comparisons of experimental 
data and hydrodynamic calculations and their interpretation, before I give 
a detailed account on recent and future developments that can be expected 
in the field. Almost every single observable of heavy ion physics is in- 
fluenced by the bulk matter, as the bulk provides the background of the 
evolution through which even the rare probes have to propagate. Now, with 
hydrodynamics, we have a tool to calculate the dynamics of this evolution. 

2. Indications of strong transverse flow 

If the microscopic constituents of a flowing medium share a common aver- 
age flow velocity v± , heavier particles in this medium gain larger momenta 
than the lighter constituents. In the case of the late stages of heavy ion 
collisions, the flowing medium is the hadronic soup consisting of hadrons 
of widely varying masses. From experimental evidence we find that the 
hadronic abundancies reflect a composition of a hadronic gas at a temper- 
ature Tchom ~ I75McV. 6 Refinements of the simple blast wave model 7,8 
assume the knowledge of a certain flow profile at the break up stage of 
the reaction as well as the shape of this surface of decoupling from which 
particles supposedly do not interact any further. In this way one achieves 
a handy formalism, 9 that can easily be used to fit to the vast amount of 
data currently available from RHIC experiments, in particular to single 
particles spectra and eventually their anisotropics. 10 Although the temper- 
ature parameter of these fits differs widely from 175 MeV 9 to 100 MeV, 10 
the common finding of these model calculation is that a vast amount of 
transverse flow (with an average velocity greater than c/2) is achieved on 
a rather short time scale of about 10 fm/c. While the stunningly good de- 
scription hints that there possibly is some physical truth underlying these 
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parameterizations, these simple descriptions of the final state do not offer 
any insight how and when such strong transverse flow develops during the 
course of the evolution. For this, we have to turn to a fully dynamical 
description of the transverse expansion stage of the reaction. 

3. Hydrodynamic formalism and phenomenology 

There are two main philosophies to model the dynamical evolution of a 
system. If the system is small and the scattering processes can be treated 
individually, one can adapt a microscopic viewpoint, treating all the colli- 
sions of the constituents individually. 11 In a very dense, strongly interact- 
ing medium, however, this approach becomes quickly impractical, partic- 
ularly if one also were to consider low-momentum transfer processes con- 
sistently. Alternatively, in this limit of high density and strong rescatter- 
ing, one can give up the particles' individual personality, and characterize 
the system in terms of density fields and continuity equations. Ideally, if 
scattering is strong enough to allow for local thermodynamic equilibrium, 
the equations of energy momentum conservation can be formulated by the 
thermodynamic equivalent, the conservation of the thermodynamic energy- 
momentum tensor, [(e + 'p)u iL u v — p g^ v ] = . The energy density e and 
the pressure p are related by the thermodynamic equation of state of nuclear 
matter, is the collective four velocity, whose time evolution we can study 
within this formalism. Once supplied with appropriate initial conditions, 
which are often taken from geometric considerations 12 or more fundamental 
models such as the saturation model 12 or the color glass condensate, 13 the 
system's evolution is fully determined. At this point, the computer takes 
over to solve for the time evolution of the thermodynamic fields. 

Fig. 1 shows the temporal evolution of the entropy density on a double 
logarithmic scale for three different positions in the fireball of a central col- 
lision, 0, 3 and 5 fm away from the fireball center. 14 The tangents drawn to 
the curve stress the occurring transition from the initial one-dimensional, 
longitudinal expansion to a fully 3-dimcnsional expansion at late stages. 
The right plot of Fig. 1 shows the radial flow profile at different times 
throughout the evolution. Striking here is the fast transition to a linear pro- 
file, which persists throughout time. Still, as matter is continually moved 
to larger radii, the mean transverse velocity of the medium increases and 
transverse flow thus continues to increase. The main part of the fireball 
matter is know to freeze-out over a rather short time frame as the rescat- 
tering rate drops sharply as a function of temperature, and therefore even 
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more so as function of time. This fact, together with the rapid development 
of a flat flow profile is at the heart of the applicability of the simple and 
the improved blast wave parameterizations. 




Figure 1. Left panel: Time evolution of the entropy density at three different points in 
the fireball (0, 3, and 5 fm from the center). Dashed lines indicate the expectations for 
pure one-dimensional and three-dimensional dilution, respectively. 14 Right panel: The 
radial velocity profile at different times during the evolution. 14 Shown is also the linear 
profile v(r) = 0.07 r/fm. 



In contrast to the overall radial flow which increases throughout the life- 
time of the system, anisotropics in the transverse flow profile are generated 
during the earliest instances of the reaction. Such anisotropics arise from 
the eccentricity in coordinate space of non-central collisions. Larger pres- 
sure gradients in the short direction (the direction of the impact parameter) 
lead to a larger transport of matter in this direction, thereby reducing the 
eccentricity and undermining the source of its own origin. 15 In the course of 
this process more particles are transported into the direction of the impact 
parameter and the mean transverse momentum of the particles is greater, 
both of which manifests itself experimentally in anisotropics v n of the trans- 
verse momentum spectra pTd ^ dydv = 27r Jd PT d y (* + En v n cos(mp)) . 
It has been shown in many studies, both microscopically 16 as well as 
macroscopically 2 ' 14 that those momentum anisotropies are generated during 
the first 5 fm/c of the reaction, where most of the matter is still at temper- 
atures exceeding the critical temperature of nuclear matter. Anisotropies 
in the transverse momentum spectra of hadrons thus originate from the 
partonic stage of the fireball and signals their collective behavior. 2 
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4. Application to experimental observations 

As mentioned before, the transverse momentum distribution at moderate 
Pt are the natural choice of observables to study collective dynamical effects 
in the medium. We use the transverse momentum spectra of identified 
particles to determine the parameters of the calculation for RHIC at 200 
GeV center of mass energy per nucleon. Those spectra are well reproduced 
applying an equilibration time of T equ = 0.6fm/c and a central fireball 
temperature of 360 MeV. 17 Although the particle ratios reflect chemical 
equilibrium at a temperature of about 175 MeV, 6 the slopes in particular of 
heavier baryons require further transverse accelertion down to temperatures 
of about 100 MeV. 18 To keep the particle ratios fixed to their value at 175 
MeV, chemical potentials seem to be dynamically generated throughout the 
evolution. 19 ' 20 ' 21 




Figure 2. Left panel: Transverse momentum spectra of Pions, Kaons and Antiprotons, 
from all four RHIC experiments combined 22 in comparison to results from hydrody- 
namics at different freeze-out temperatures (and different initilizations, see 17 ). Right 
panel: Momentum anisotropics of Pions, Kaons and Lambdas. Results from the STAR 
experiment 23 compared to results from hydrodynamics. 



Although the good reproduction of the spectra of central collisions is 
reassuring, it has been achieved to a great deal by adjusting the parameters 
of the calculation. Important predictions come into play in non-central col- 
lisions. Changing centrality, which introduces a breaking of the azimuthal 
symmetry inherent in central collisions, renders essentially all observables 
sensitive on the azimuthal angle of observation. This literally opens up an 
entirely new dimension of observable space. However, for fully dynamical 
simulations no additional parameters are introduced when stepping out to 
explore observables in this new dimension. The modified initialization for 
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non-central collisions is stricktly determined by the geometry of the under- 
lying initialization scenario. The only appearing parameter - the impact 
parameter b - is estimated experimentally by giving the range of centrality 
for certain collisions. 

We already mentioned the great significance of momentum anisotropies 
as they provide a signal from the first few fm/c of the collision. Microscopic 
studies 11 in comparison with experimental data show that the fireball con- 
stituents must undergo an incredible amount of rescattering among them- 
selves. Hydrodynamic calculations, which exploit the limit of negligible 
mean free pathlcngths, provide an impressive account of a large collec- 
tion of momentum anisotropies. 18 They give a good overall description of 
momentum anisotropies up to transverse momenta of 2 GeV (Pions) and 
even 4 GeV for Baryons. The centrality dependence of the hydrodynamic 
elliptic anisotropy compares well with experimental data as long as the im- 
pact parameter does not get too large. Finally, the mass characteristics, a 
flatter onset of momentum anisotropies at small transverse momenta, 25 is 
confirmed by experiment in stunningly good agreement. A recent compi- 
lation of momentum anisotropies of heavy particles is reprinted in Fig. 2 
(right panel). 23 Still more recently the momentum anisotropies of Cascades 
and Omegas have become available. 24 Again, those anisotropies have been 
stunningly large, as predicted 3 years ago by hydrodynamic calculations. 25 
As the multistrange resonances are believe to not interact strongly in a 
hadron gas, this is the clearest signal that strangeness anisotropy is gen- 
erated in the partonic stage of the reaction. It appears that all quark fla- 
vors share the same flow anisotropy on the partonic level. This fact could 
be most prominently (dis-)proved by investigating momentum anisotropies 
of higher order and their ratios. 26 That the momentum anisotropies are 
as large as observed is a strong indication that the thermalization of the 
medium is achieved very rapidly. 27 The fact that hydrodynamics overesti- 
mates anisotropies at large transverse momenta and in peripheral collisions 
is well understood, considering that high px particles escape the fireball too 
rapidly to follow the collective motion of the bulk, and that small systems 
do not provide enough rescattering in the limited volume. Here, it is also 
important to remark that while hydrodynamics does not seem to be fully 
applicable at the lower beam energies of the SPS in peripheral collisions, 
there might still be a fair chance for its quantitative application in central 
collisions. This could be in particular true for the early, plasma part of the 
reaction, whereas the later hadronic stage might drop out of local thermal 
equilibrium. Only detailed quantitative comparisons of hydrodynamic cal- 
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culations and the most recent results 29 from SPS experiments can answer 
this question and address important topics close to thermal equilibrium and 
just above the transition temperature. 

5. Current trends and future requirements 

Most of the hydrodynamic calculations that are applied to experimental 
data at present include some great simplifications in order to keep the 
numerical efforts to a reasonable degree. One of these is the strong bias 
to investigations of observables at midrapidity. Although the collectivity 
observed at midrapidity is clearly produced during the evolution of the 
produced fireball and therefore contains clean probes of the dynamics, a 
broader view of the full collision volume is highly desirable, in particular in 
the light of recent data extending the analysis of momentum anisotropics to 
large rapidities. Fully three-dimensional calculations have been successfully 
performed in recent years. 30 However not too much is gained by applying 
the ideal hydrodynamic calculations at large forward rapidities as the sys- 
tem is thought reach thermal equilibrium only at later times 31 or not at 
all 32 in regions far away from the center of mass of the reaction. 

Nevertheless, there is a lot to be gained from investigations of hydrody- 
namic calculations at non-zero rapidity. Due to the (net)-baryon content 
of the initial nuclei, one can expect that the (net) baryon-number density 
increases with increasing rapidity. As one leaves the 'clean' midrapidity 
region with its antibaryon to baryon ratio of about 0.75, one gets closer to 
the initial nuclei region where the antibaryon content effectively vanishes. 
Thus, at larger rapidity, a hydrodynamical system evolves according to the 
equation of state at larger baryon chemical potential than at mid-rapidity. 
This offers the exciting possibility of tracing signals of the expected tricrit- 
ical point at finite rapidity. The effect of the tricritical point on the hydro- 
dynamic evolution of RH1C collisions is currently explored by Nonaka. 33 ft 
was found that such a point acts as an attractor for the thermodynamic 
paths of a constant fraction of entropy and baryon-density. Clearly, due 
to the large fluctuations that appear around the critical point, the ideal 
hydrodynamic treatment will become invalid, but it will be interesting to 
explore how the system drops out of equilibrium, and which quantitative 
experimental measurements will be predicted. 

The less efficient energy density production and deposition at large ra- 
pidity eventually leads to a breakdown of the assumption of full thermal- 
ization and an ideal expansion. Viscous effects will become more and more 
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important in smaller and more dilute systems. It is of fundamental interest 
to get estimates of the viscosity of the QCD plasma state. Viscous hydrody- 
namics can deliver valuable information in this direction. The realization 
of fully dynamical hydrodynamical calculations is however difficult. The 
viscous terms have to be solved dynamically, increasing the number of dif- 
ferential equations by up to 5 (depending on the symmetries one exploits 
explicitly). Although difficult to solve, 34 first results on the dynamical evo- 
lution of central collisions have recently been presented. 35 The message of 
these dynamical studies is that viscous matter sticks together for a longer 
period, just to blow up quicker in the late stages. The left panel of Fig. 
3 shows plots of the energy density as a function of the transverse coor- 
dinate for a viscous calculation compared to the ideal Euler calculation. 35 
Although the densities change slower at first, velocities build up faster due 
to a restructuring of the pressure components, which is also observed in the 
case of only transverse thermalization and longitudinal free streaming. 36 Fi- 
nally, in the late stage of the reaction, the left over matter dissipates more 
rapidly, due to the larger velocities that were achieved. The net result is 
that the system stays longer in the hot phase of its reaction and disperses 
faster through the late stages. The late stages, for which viscosity become 
increasingly important, thus have a shorter lifetime and freeze-out occurs 
more rapidly. Comforting is the observation that viscosity does not lead to 
self-enhancing features! 




Figure 3. Left panel: Radial change of the energy density at various times in the ideal 
calculation and in a viscous calculation. 35 Right panel: Changes on the final state mo- 
mentum anisotropy due to deviations from the ideal of momentum space distributions. 37 
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In this context it is also interesting to cast doubt on the assumption 
of having particles distributed in momentum space according to the simple 
ideal distribution laws of Fermions and Bosons. It has been shown 37 that a 
large deviation from these distributions smears out the directed collective 
motion observed in the momentum anisotropies and is thus not supported 
by the data (see right panel of Fig. 3.) 

6. Summary 

We have learned a lot about the collective phenomena exhibited in RHIC 
data of the first few years. The data clearly show many characteristics 
of rapid thcrmalization and an extended stage of hydrodynamic expansion 
of the fireball into the surrounding vacuum. On a linear timescale, the 
hydrodynamic stage is the most significant stage of the reaction, thereby 
influencing many other relevant observables. Hydrodynamics is a clean 
(meaning well defined) tool, with only a small number of parameters and 
easily reproducible by anybody having access to a computer. The two most 
important facts learned from hydrodynamics and the experimental data so 
far is the need for rapid thermalization (in a time of less than 1 fm/c) 
where matter obtains a large transverse push due to a hard equation of 
state (p ~ e/3). Many quantitative questions are still unanswered: What 
is the frame of error around the equilibration time, what different kinds of 
equations of state and phase transitions would still account for the data, 
how does one feed back on the other? What is the role of viscosity in the 
expansion? Why does viscosity of high temperature QCD matter seem to 
be so low as to allow for an ideal hydrodynamic description of the data? 
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